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ABSTRACT 

Context. Two incorrect arguments against MOND in elliptical galaxies could be that the equivalent circular velocity curves tend to 
become flat at much larger accelerations than in spiral galaxies, and that the Newtonian dark matter halos are more concentrated than 
in spirals. 

Aims. Here, we compare published scaling relations for the dark halos of elliptical galaxies to the scaling relations expected for 
MONDian phantom halos. 

Mettiods. We represent the baryonic content of galaxies by spherical profiles, and their corresponding MONDian phantom halos by 
logarithmic halos. We then derive the surface densities, central densities, and phase space densities and compare them with published 
scaling relations. 

Results. We conclude that it is possible to get flat circular velocity curves at high acceleration in MOND, and that this happens for 
baryonic distributions described by Jaffe profiles in the region where the circular velocity curve is flat. Moreover, the scaling relations 
of dark halos of ellipticals are remarkably similar to the scaling relations of phantom halos of MOND. 
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1. Introduction 

While data on large scale structures point towards a Universe 
dominated by dark matter and dark energy, e.g. |Komatsu et al. 
(|201 1 ), the nature of these is still a deep mystery (e.g., Frieman 
[etal.,2008, Wiltshh-e 2008, Bertone,2010,.Ki-oupa et al..2010| l. In 
this context, it is good to keep in mind that this conclusion essen- 
tially relies on the assumption that gravity is correctly described 
by Einstein's General Relativity in the extreme weak-field limit, 
a regime where the need for dark matter itself prevents the the- 
ory from being tested. Until this double dark mystery is solved, 
it is thus worth investigating alternative paradigms and their im- 
plications. 

For instance. Modified Newtonian dynamics ( |Milgrom|1983| 
MOND) naturally explains various spiral galaxy scahng rela- 
tions ( |Tully & Fis"herl[T977| [McGaugh efaL] [20001 [McGaugh 
|2004| l. The existence of a very tight baryonic Tully-Fisher rela- 
tion for disk galaxies (McGaugh 2005 ; Tra chternach et al.|2009| l 
is for instance one of the remarkable predictions of MOND. The 
coiTesponding relation for early-type galaxies is much more dif- 
ficult to investigate because they are pressure-supported systems, 
and the equivalent circular velocity curves determined from the 
velocity dispersion profiles suffer from the well-known degen- 
eracy with anisotropy. However, some studies circumvented this 
problem: for instance, [Kronawitter et al" ( 2000[ l used data on 
21 elliptical galaxies to construct non-parametric models from 
which circular velocity curves, radial profiles of mass-to-light ra- 
tio, and anisotropy profiles as well as high-order moments could 
be computed. This led [Gerhard et al.| ( [200T] hereafter GOl) to 



publish benchmark scaling relations for ellipticals. It was e.g. 
shown for the first time that circular velocity curves tend to 
become flat at much larger accelerations than in spiral galax- 
ies. This would seem to contradict the MOND prescription, for 
which flat circular velocities typically occur well below the ac- 
celeration threshold gq ~ 10"^ cm s"^, but not at accelerations of 



the order of a few times ao as in ellipticals. Also Thomas et al 



(2009 hereafter T09) published scaling relations for dark matter 



halos of 18 Coma galaxies, using similar prescriptions as GOl. 
We remark that GOl employed spherical models while the mod- 
els of T09 are axisymmetric. 

Not many studies have considered the predictions of MOND 
in elliptical galaxies. Milgrom (1984) showed that pressure- 
supported isothermal systems have finite mass in MOND with 
the density at large radii falling approximately as r""*. It was also 
shown that there exists a mass-velocity dispersion relation of the 
form (M/10"Mq) » (0-^/100 kms"')'* which is similar to the 
observed Faber-Jackson relation (Sanders 2000 2010), and that, 
in order to match the fundamental plane, MOND models must 
deviate from being strictly isothermal and isotropic: a radial orbit 
anisotropy in the outer regions is needed ( Sanders'2000': 'Cardone 
et al. 20lT] ). [Tiret et al.[ ([2007l) and [ Angus eTar., ( 2008 ) also an- 
alyzed the distribution of velocity dispersion of PNe on scales 
of 20 kpc, and of satellites on very large scales of the order of 
400 kpc around red isolated ellipticals, showing that MOND al- 
lowed to fit both scales successfully. 

Hereafter, we make general remarks on the properties of 
spherical galaxies within MOND, and their scaling relations. We 
first point out a remarkable property of elliptical galaxies ex- 
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hibiting a flattening of their circular velocity curve at small radii: 
such a flattening in the intermediate gravity regime is actually 
generated by a baryonic density distribution following a JafFe 
profile in these parts of the galaxies. We then further show that 
the observational scaling relations for the dark halos of the ellip- 
tical galaxy sample by GOl are strikingly similar to the theoreti- 
cal "phantom" halos of MOND (i.e. the halo that would produce 
in Newtonian gravity the same additional gravity as MOND), 
with one remarkable exception: MOND predicts that the product 
of the central density with the core radius should be constant, as 
recently observed for spiral galaxies ( Donato et al. 2009 Gentile 
let al.|20091 l. 

2. Flat circular velocity curves and the Jaffe profile 



Although it has been argued that some ellipticals do not need 
dark matter or enhancement of gravity ( [Romanowsky 
there are many counter-example s ([Magorrian & 
200T1 IRichtler et al. 2004; Schuberth et al.||2006 



any 



et al. 2003 



Ballantyne| 



Kumar et al. 2007). Such a recent example is the elliptical galaxy 
NGC 2974 where the presence of an HI disk allowed a more or 
less direct measurement of circular velocities ( [Weijmans et al. 
2008). There is also evidence that elliptical galaxies exhibit flat 



circular velocity curves , but that, contrary to spiral galaxies, this 
happens in the inner regions where g > (e.g., GOl, Weijmans 
et al. 2008). Such a flattening of circular velocities is a priori not 
expected in the strong to intermediate gravity regime in MOND, 
and poses the question of how to analytically interpret it. 

In the intermediate gravity regime, the transition from 
Newtonian to MONDian dynamics is described by the fi- 
function of MOND. Many concordant studies have recently 
shown that, in spiral galaxies, the "simple" transition of Famaey 



& Binney ( 2005]l is a good representation of the data (Gentile, 



Famaey, & de Blok| 2011 for an extensive discussion). In a 
spherical system, with this simple transition, the enclosed (bary- 
onic) mass Mjyiir) needed to produce the same gravitational 
potential in MOND as the (baryonic-i-dark) mass M^ir) in 
Newtonian gravity is: 



\MN{r) r^aoj 



(1) 



In a region where the circular velocity is constant vv = V (even 
if g > ao), one can write M^ir) - V^r/G, and thus after some 
algebra 



y4 

MM(r) = — ■ 



fl()G r + V^/ao 



(2) 



Remarkably, this enclosed ma ss profile corresponds precisely 
to a JafFe profile (JafFe 1983 1 with scale-radius rj - V^/ao 
(meaning that the acceleration is ao at rj), and with total mass 
Mtot = y"*/(fl()G). Indeed, as the enclosed mass Mm(/) - 
A^M('"o)+47r r p(R)R^dR, this enclosed mass profile corresponds 
locally to the density profile: 



P(r) 



47r 



(r + rjY 



(3) 



Mtot/rj = ao/G. This profile is of course not valid for the very 
inner parts of an elliptical galaxy, where V is not constant. Let 
us also note that (i) it was already known that a Jaffe profile 
produces a flat circular velocity curve at r rj in Newtonian 



gravity, which MOND generalizes to radii r ~ ry, (ii) Mtot does 
not necessarily have to be the real total mass of the galaxy, as 
the JafFe profile fit to the density distribution could have a cut- 
off in the outer parts. In that case, the constant circular velocity V 
would actually Fall slightly above the prediction From the bary- 
onic Tully-Fisher relation oF spiral galaxies. Interestingly, this is 
precisely what is observed For the GOl sample oF ellipticals. 

The Fact that elliptical galaxies can exhibit (equivalent) cir- 
cular velocity curves that are Flat in the intermediate gravity 
regime is thus analytically understood in MOND by the Fact 
that the outer regions oF ellipticals can be approximated by a 
JafFe profile with a large scale-radius, i.e. in regions well within 
the intermediate gravity regime rather than in the deep-MOND 
regime. These flat circular velocity curves would have been im- 
possible with exponential density profiles (as encountered in spi- 
ral galaxies), meaning that the Fact that circular velocity curves 
become flat quicker in ellipticals does not come as a surprise in 
the context oF MOND. 

This finding looks like an interesting possibility to devise 
new tests oF MOND based on photometry. However, in reality 
it might be difficult: not many spherical galaxies with a pre- 
cisely measured density profile are dynamically investigated out 
to large radii, and have enough tracers to measure the higher 
order moments and constrain the anisotropy. Moreover, light 
might not trace the baryonic mass precisely. As an example, 
the circular velocity in NGC 2974, which can be traced by an 
HI disk, becomes constant at around 5 kpc and has the value 
300 kms"', which would correspond to a Jaffe scale radius of 
23 kpc. Unfortunately, NGC 2974 is neither spherical nor does 
its photometry reach large radii so that it does not serve well as 
a test object. In any case, Weijmans et al. (2008, their Fig. 20) 
showed that the reverse procedure (going from the density to the 
circular velocity curve) leads to a very good fit. 

3. Dark matter scaling relations for phantom halos 
of ellipticals 

We now apply the reverse procedure, and check whether the 
phantom halos predicted by the simple transition of MOND 
(Famaey & Binney 2005) comply with the observational scal- 
ing relations of dark halos of ellipticals. As stated above, Jaffe 
profiles are not good descriptions of the very inner parts of ellip- 



ticals. We hereafter rather choose Hernquist profiles (Hernquist 



1990) to represent the baryonic content of ellipticals: these are 
realistic enough and allow for an exhaustive exploration of their 
properties without varying too many free parameters. Such a 
Hernquist-model is described by its total mass M and scale- 
radius rH- The profile of the Newtonian circular velocity curve 
then reads 



VNir) 



GMr 



\ (r + rn)^ 



(4) 



where the scale-radius rn of the Hernquist-model is related to 
the effective (half-light) radius by R^jf - 1.815 rn- 

Adopting the simple transition formula between the 
Newtonian and the MONDian regime, one finds for the MOND 
circular velocity 



with the characteristic surface density (see also Milgrom 1984) vm 



^Jvl{r)/2 + ^v^(r)/4 + vl(r)aor, 
DNDii 

i-'phantom — sj^' 



and the MONDian phantom halo has the circular velocity 



vUr)/4 + vUr)aor-viir)/2 



(5) 



(6) 
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Table 1. The table shows for baryonic Hernquist profiles with mass and effective radius described by the first two columns, the 
corresponding parameters of the MONDian "phantom" halo represented by a logarithmic potential fitted from the center to two 
effective radii of the baryonic profile. The columns are the baryonic mass, the luminous effective radius, the core radius ro and 
the asymptotic velocity vq of the log-halo, its surface density, central density and phase space density, the latter as defined by 
GOl.The last column column gives the acceleration in units of ao at a radius of 2 Re/j- for each Hernquist model. The predictions 
for MONDian halos in the previous columns are valid only for galaxies embedded in an external field smaller than this value. 



baryonic mass [Mg] 


Reff [kpc] 


ro [kpc] 


vo[kms '] 


SlMJpc^] 






acc.[ao] 


10'^ 


14.1 


8.83 


244 


374 


0.04 


8.82X 10-'-' 


1.53 


8x 10" 


11.8 


7.60 


228 


379 


0.05 


1.19 X lO-** 


1.55 


5x 10" 


8.06 


5.20 


193 


393 


0.076 


3.02 X lO^** 


1.91 


2x 10" 


3.84 


2.87 


146 


411 


0.143 


1.3 X 10^' 


2.96 


10" 


1.47 


1.24 


99 


438 


0.35 


1.02 X lO^*^ 


8.44 


5 x 10'" 


1.25 


1.04 


90 


431 


0.41 


1.6 X 10^* 


6.08 



To enable the comparison with the scaling relations of GOl, 
where the dark matter halos are adopted as logarithmic halos, we 
fit Vphantom('") to the circular velocity viog(r) of a logarithmic halo 
with asymptotic circular velocity vq and core radius ro: 



Viog('") = vor/ ^ 



(7) 



The fits are performed within the inner two effective radii The 
fitted central dark matter density is then given by 



Po = 3(yo/ro)V(4;rG). 



(8) 



The characteristic central phase space density is defined (see 
GOl) as 

fps = ^"^Polvl- 

The characteristic surface density within ro is then also defined 
as (see also Donato et al. 2009): 



(10) 



Table [T] lists these fitted parameters for six baryonic Hernquist 
masses over a large mass range. The combinations of the masses 
and effective radii in Table [T] follow equation (5) of GOl (in ac- 
cordance with the fundamental plane), where we transformed 
their luminosities into masses by using M/Lg-Bi for all galaxy 
baryonic masses. 

Fig[T]shows the values of the fitted dark halo parameters de- 
rived from applying MOND to the baryonic Hemquist-profiles, 
together with the observational scaling relations given by GOl 
(dotted lines) and T09 (dashed lines). The upper panel shows 
the characteristic phase space density, the middle panel the cen- 
tral volume density, and the lower panel the characteristic sur- 
face density. Let us note that the plotted relations are indica- 
tive only, since the data (Fig. 18 of GOl and Figs. 1 and 4 of 
T09) show a very large scatter even when logarithmically dis- 
played. However, within this observational uncertainty, it is re- 
markable that some features are perfectly reproduced, particu- 
larly the slopes of the phase space density and of the central 
volume density as a function of baryonic mass (given the ob- 
servational scatter, the almost perfect reproduction of the central 
volume density of GOl might of course be partly coincidental). 

As first emphasized by GOl, the phase-space density values 
are at a given mass higher than in spirals, which means that un- 
der the A Cold Dark Matter paradigm, dark halos of ellipticals 



' Let us note that these fits are not particularly good: the circular 
velocity curve V(.(r) = vor/(ro + r) would have provided better fits, but 
the core radius of the corresponding halo would then be systematically 
smaller with respect to GOl 



cannot be the result of collisionless mergers of present-day spi- 
rals, but must have been assembled at a very early time, when 
the cosmological density was higher In MOND this is of course 
not necessarily the case, as the phase-space argument does not 
apply to phantom halos. 



One also notes a remarkable exception to the scaling rela- 
tions: the fitted characteristic dark matter surface density S is 
fully independent from the Hernquist parameters, and it is sys- 
tematically lower than in GOl and T09. We emphasize that this 
constancy is not related to the special relation of mass and effec- 
tive radius. Varying R^j f by a factor of two at a given mass does 
not change the constant surface density significantly. This pre- 
diction of MOND thus brings the value closer to the (also con- 
^gy Slant) value of S observed in spiral galaxies, log5 =2.1 ( [Donato 



et al. 2009). Let us note that MOND also predicts the observed 
constant value of S in spirals, which is somewhat lower because 



(i) spirals are a bit deeper into the MOND regime (Milgrom 
12009) and (ii) their flattened baryonic profiles lead to a some- 
what higher Newtonian gravity at a given mass, and in turn a 
somewhat lower MOND contribution to the phantom halo. 

On the first glance one might interpret this constancy and 
the other scaling relations as a clear signature of MOND in el- 
lipticals: however, CDM may also predict that the surface den- 
sity within the scale radius of NFW halos weakly depends on 
dark matter total mass (j Boyarsky et al.|2010[ l. For spiral galax- 
ies, this is of little interest as it is known that cuspy profiles of- 
ten do not fit rotation curves ( Famaey & Binney 2005 ; de Blc^ 
2010 Gentile et al.|2005" i, the mystery then being how to erase 



the cusp by the feedback from the baryons while keeping the 
product po'"o constant. In elliptical galaxies, the situation is less 
clear as NFW profiles often do fit the data equally well as cored 
profiles ( |Schuberth et al.,2010) . We thus fitted NFW profiles to 
the same MONDian phantom halos and found a perfect agree- 
ment. The question remains whether these NFW-halos are "cos- 
mological" or in other words, fulfill the relation between virial 
mass and concentration predicted by cosmological simulations. 
Fig|2] displays for our Hernquist masses the resulting concen- 
trations of the NFW-halos (open circles) corresponding to the 
MONDian phantom halos, while the triangles show the concen- 
tration values expected from the equation (9) of |Macci6 et al. 



(2008), using 200 times the critical density as the mean density 
within the vmal radius (standard cosmology: h=0.7, Q,„ - 0.3, 
Qa = 0.7). One concludes that for high masses the MONDian 
phantom halos are not distinguishable from cosmological NFW 
halos, given also that the simulations predict considerable scat- 
ter. For smaller masses the difference between MONDian phan- 
tom halos and NFW cosmological halos is larger. 
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Fig. 1. The figure shows the surface density, central density, and 
the central phase space density logarithmically (see the text for 
more explanations) of the phantom dark halos (circles) for dif- 
ferent baryonic Hernquist masses from Table[T]together with the 
relations given by GOl (dotted lines) and T09 (dashed lines). 
Note that these parameters observationally exhibit a very large 
scatter around the mean relations. Note also that the relations of 
T09 are not given explicitly in their paper but have been con- 
structed from their Table 3 omitting galaxies with young stellar 
cores. 



4. External field effect 

Due to the non-linearity of MOND and its associated breaking 
of the Strong Equivalence Principle, a MONDian stellar system 
embedded in an external gravitational field (EF) stronger than 
its own internal field behaves in a quasi-Newtonian way, with an 
effectively higher gravitational constant (Milgrom 1983[ Famaey 
et al.|2007) . Most of the sample galaxies are located in clusters 
or groups where the EF might have an influence. |Wu et al^ ( |2010| l 
for instance showed how the EF can lead to the lopsidedness of 
an originally axisymmetric non-isolated galaxy. 

While it is beyond the scope of this research note to evaluate 
in detail the EF in the present sample, a very rough estimation 
is presented in Fig. |3] which plots for the Virgo and the Coma 
cluster the accelerations based on the extrapolations of the mass 
models cited in the figure caption (these extrapolations are only 
meant to give an order of magnitude estimate, but should not be 




11 12 
log mass Hernquist 

Fig. 2. The figure shows for our 6 Hernquist masses the con- 
centration parameters of the associated NFW-halos, if the 
MONDian phantom halos are fitted by NFW profiles (open cir- 
cles). The triangles are the concentration parameters expected 
from the relation quoted by Maccio et al. ( 2008| l 
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Fig. 3. This plot estimates the external field acting on Virgo 
(dashed) and Coma (dotted) galaxies. Abscissa is the pro- 
jected distance in Mpc from M87 and NGC 4874, respectively. 
Ordinate is the acceleration in units of ao. Adopted distances 
for Virgo and Coma are 15 Mpc and 100 Mpc, respectively. 
The values for Virgo are generated by an extrapolation of the 
mass model for M87 of McLaughlin ( 1999 ). Some GOl galaxies 
are indicated by their NGC numbers. The values for Coma are 



generated by using the NFW dark halo from Lokas & Mamon 
( |2003| l. Small open circles are the 18 T09 galaxies whose pro- 
jected distances are taken from Godwin et al. ( 1983, ). The com- 
parison with Table [T] shows that the EF is small. 



taken as rigorous models). This can be compared with the in- 
ternal accelerations at 2 Re/f for the Hemquist models in TabfT] 
Indicated are the projected distances of galaxies in the Virgo and 
Coma region. The positions of the Virgo galaxies correspond to 
the middle points of their NGC numbers, while the Coma galax- 
ies are plotted as small open circles. One concludes that the EF 
should have no influence in the two samples at the galactocentric 
distances which we consider. 
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5. Conclusion 

Here we showed that (i) in MOND, galaxies exhibit a flatten- 
ing of their circular velocity curve at high gravities (g > ao) 
if they are described by a JafFe profile with characteristic sur- 
face density oq/G in the region where the circular velocity is 
constant (since this is not possible for exponential profiles, it is 
remarkable that such flattenings of circular velocity curves at 
high accelerations are only observed in eUiptical galaxies); (ii) 
the phantom halos of ellipticals predicted by MOND (i.e., the 
dark halos that would produce in Newtonian gravity the same 
additional gravity as MOND) can be fitted by logarithmic halos 
which perfectly reproduce the observed scaling relations of el- 
lipticals for phase-space densities and central volume densities 
po; (iii) these halos have a constant characteristic surface den- 
sity poro; (iv) contrary to spirals (for which there are more data 
in the very central parts), the phantom halos of ellipticals can as 
well be fitted by cuspy NFW halos, the concentration of which is 
in accordance with the theoretical predictions of ACDM for the 
highest masses, but in slight disagreement for baryonic masses 
smaller than lO^M©: a modern, large, sample of elliptical galax- 
ies, which are dynamically well investigated out to large radii 
and cover a large range of masses, will thus be required to get 
discriminating power But in any case, and whatever the true 
physical reason for it, it is remarkable that a recipe (MOND) 
known to fit rotation curves of spiral galaxies with remarkable 
accuracy also apparently predicts the observed distribution of 
"dark matter" in elliptical galaxies. 
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